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ABSTRACT 

Type la supernovae (SNe la) can be used as the standard candlelight to determine the cosmological 
distances because they are thought to have a uniform fuel amount. Recent observations of several 
overluminous SNe la suggest that the white dwarf masses at supernova explosion may significantly 
exceed the canonical Chandrasekhar mass limit. These massive white dwarfs may be supported by 
rapid differential rotation. Based on single degenerate model and the assumption that the white 
dwarf would differentially rotate when the accretion rate M > 3 x 10 Mgyr -1 , we have calculated 
the evolutions of close binaries consisting of a white dwarf and a normal companion. To include the 
effect of rotation, we introduce an effective mass M e g for white dwarfs. For the donor stars with two 
different metallicities Z = 0.02 and 0.001, we present the distribution of the initial donor star masses 
and the orbital periods of the progenitors of super-Chandrasekhar mass SNe la. The calculation results 
indicate that, for an initial massive white dwarf of 1.2M©, a considerable fraction of SNe la may result 
from super-Chandresekhar mass white dwarfs, but very massive (> 1.7M Q ) white dwarfs are difficult 
to form, and none of them could be found in old populations. However, super-Chandrasekhar mass 
SNe la are very rare when the initial mass of white dwarfs is 1.0M©. Additionally, SNe la in low 
metallicity environment are more likely to be homogeneous. 

Subject headings: binaries: close — stars: evolution — supernovae: general — white dwarfs 



1. INTRODUCTION 

Type la supernovae (SNe la) are generally believed 
to be thermonuclear explosions of accret ing carbon- 
oxygen white dwarfs (|Hovle fe Fowler lll960| ) when their 
masses reach the cha ndrasekhar mass of ~ 1.4M 
( Chandrasekh arl I1931D . Due to this uniform progeni- 
tor mass ( Mazzzli et al. I I2007T ) . most SNe la present a 
good correlation between the peak brigh tness and the 
width of the light curve (Phillips 1993), so it is pos- 
sible to use them as the stand ard candlelight to deter- 
mine the cosmological d istances (jRiess et al. H9Ml200l 
iPerlmutter et al. II 19991) . 

There exist double degenerate model 



1981 



a nd single 

(|Whelan fc Iben I 11973 : iNomoto I 



Ijlben fc Tutukov I [l98l fW ebbink 
degen erate model 

19821) for the proge nitors of SNe la (for a review see 
Branch et al. 1 119951 ) . Based on the single degenerate 
model, the evolutions of binary systems consisting of 
an accreting white dwarf have been widely explored 
by many authors (lHachisu. Kato fc Nomotol 119961: 
Li fc van den Heuvell 119971: lHachisu. Kato fc Nomotol 

19991 
2004; 



also suggested white dwarfs with super-Chandrasekhar 
masses as the progenitors of some rare nickel-rich SNe 
la (MNi > 0.8M©). These ma ssive white dwarfs ma y 
be supported by rapid rotation (jYoon fe Langer I f2005^ . 
or origin from the merger of two massive w hite dwarfs 
(|Tutukov fe Yungelsonlli99l iHowelfl l2001h . More re- 
cently, two other possibly overluminous SNe la (SN 
2j^z. lfflck"en et al. I ((20071) and SN 2007if. lYian"et al. I 
(2007)) were reported. 

It is interesting to explore the properties of the pro- 
genitors of overluminous SNe la like SNLS-03D3bb in 
the single degenerate model. Recently, some authors 
have taken account of the influence of rotation on the ac- 
creting white dwarf in t heir works (lUenishi et al.l 120031 : 
Saio fc Nomoto I l2004t lYoon fc Langer I I2004L l2005f >. 
Yoon fc Langer I (|2004l ) found that a white dwarf with 
accretion rate of > 3 x 10 _7 M Q yr _1 may rotate differ- 
entially. The maximum equilibrium mass for a secularly 
stable white dwarf with diff erential rotation is given by 
dShapiro fc Teukolskv lll983h 
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The supernova SNLS-03D3bb (2003fg) was discovered 
on 24 April 2003, and observed t o be 2.2 times over - 
luminous tha n a n ormal SN la (|Astier et al. I 120061 ). 
iHowell et al~l (|2006f ) inferred the mass of 56 Ni ~ 1.3M , 
which indicates that the mass of th e white dwarf at the 
moment of explosion is ~ 2.1M© 1 . IWang et al. I (|2008l ) 



1 IHillebrandtTS im & Ropkc (2007) suggested an alternative 
model to account for the high brightness of SNLS-03D3bb, in which 
an off-center ignition of nuclear burning in a Chandrasekhar mass 



M max « 2.5( — ) 2 M Q , 
Mc 

where fj, e is the mean molecular weight p er electron. This 
valu e is in line with the results de rived bv lDurisen I (|1975l ) 
and iDurisen fc Imamura 1 ([T981) via detailed numerical 
calculations. 

Consider ing the angular momentum tr a nsfer b y mass 
accretion, lYoon. Langer fc Scheithauerl (|2004l ) simu- 
lated the evolution of helium-accreting CO white dwarfs 
for different mass accretion rate, and found that the ro- 
tation can help white dwarfs grow in mass by stabilizing 
helium shell burni ng. Based on th e rigid ly rotating pro- 
genitor model, iDommguez et al.l ((2006) computed the 

white dwarf occurs and drives a lop-sided explosion 
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evolution of a rotating CO white dwarf accreting CO-rich 
matter. Their results show that more massive progeni- 
tors result in higher 56 Ni mass and explosive luminosity, 
and more massive white dwarfs at explosion. 

Assuming that overluminous SNe la originate from bi- 
nary systems consisting of an accreting white dwarf with 
super-Chandrasekhar mass, in this paper we investigate 
the distribution of the initial donor star mass and or- 
bital period of the progenitor binaries. In section 2, we 
describe the input physics in stellar and binary evolution 
calculations. Numerically calculated results for the evo- 
lutionary sequences of white dwarf binaries are presented 
in section 3. In section 4 we summarize the results with 
a brief discussion on the limitation of this study. 

2. INPUT PHYSICS 

In this work, we study the evolution of close binaries 
consisting a CO white dwarf (of mass M wd ) and a nor- 
mal companion (of mass M d ) with population I and II 
metallicities (Z — 0.02 and 0.001, respectively). The 
white dwarf binary results from a binary consisting of 
two main-sequence stars, in which the more massive star 
evolves more rapidly, and becomes the CO white dwarf. 
With nuclear evolution, the secondary star starts to fill 
its Roche lobe and transfer hydrogen-rich material to the 
white dwarf. Our calculations begin at this stage. 

The accreted material on the white dwarf will be 
heated and compressed, finally leading to nuclear burn- 
ing. If the burning process is unstable, part of the mass 
is ejected from the white dwarf. The accumulated effi- 
ciencies for hydrogen and helium burning are denoted to 
be an and an c , respectively. Altho ugh calculations of 
an an d QH e have been done b y, e.g . jKovetz fc PrialnikI 
(|1994h and lKato fc Hachisul (|2004h . they are limited to 
non-rotating white dwarfs, and in principle, inadequate 
to determine the final mass of H accreting rotating white 
dwarfs. To explore the effect of stellar rotation on the 
thermal response of accreting CO white dwarfs, we in- 
clude the lifting effect in the hydrostatic equilibr i um ac - 
cording to the prescriptions in lDoimnguez et al~l (|1996h . 
and introduce an effective mass M e g of the white dwarf 
by taking account of the centrifugal force. We then adopt 
the corresponding values o f Qh and ctHe for this effective 
mas s from the formulae in lHan fc Po dsiadlow ski I (|2004l ) 
and iKato fc Hachisul (|2004h . The method can be de- 
scribed as follows. We divide the surface of the white 
dwarf into three zones with different ranges of the polar 
angle: the equatorial zone (9 = 60° — 120°), the middle 
zone (9 = 30° - 60° and 120° - 150°), and the polar zone 
(9 = 0° - 30° and 150° - 180°). Assume rigid body rota- 
tion and neglect any deformation of the white dwarf, the 
effective mass M e g of the white dwarf in the zone with 
polar angles between 9\ and 9i satisfies 



GM eS _ GM wd Ijl 27r^Rhm 3 9d9 
fjj* 2irR 2 sin6dd 



R? 



R 2 



(2) 



where R and ui are the radius and angular velocity of 
the white dwarf, respectively. In Eq. (2), the radial com- 
ponent of the centrifugal force is averaged by an area- 
weighted mean. 

We further assume that each zone accretes the trans- 
ferred material at a rate proportional to its area, i.e. 

its accretion fraction /; = f 9 * 2ir R 2 sm9d9 / 4n R 2 , and 



/i = 0.5, 0.366, and 0.134 for the equatorial zone, mid- 
dle zone, and polar zone, respectively. Replacing M wd 
with M c ff, we can simply obtain the H and He accumu- 
lated efficiencies («H,i and QfHe,i) for different zones on 
the surface of the white dwarf. Summarize the above 
prescriptions, the mass growth rate of the white dwarf 
can then be written as 



M wd = 



^ a H ,iaHe,i/i-Md 



(3) 



where Md is the mass transfer rate from the donor star. 

We also calculate the spin evolution of the white dwarf, 
neglecting the possible interaction between the mag- 
netic field of the white dwarf and the accretion disk. 
The rotation of the white dwar f may be related to the 
rotation of its progenito r star (He ger fc Langer 1 120001 : 
iMaeder fc Mevnet I l2000T >. but is more likely to be at- 
tained during mass transfer, as the transferred mate- 
rial from the donor star carries a large amount of an- 
gular momentum, wh i ch can cause spin-up of the white 
dwar f (iDurisen 1 119771: IRitter Ill985t iNaravan fc Pophm I 
Il989t lLanger et al. I I2000T ) . Generally the wh ite dwarf 
will spin up along w ith accretion from a disk (|Durisen I 
119771 : IRitter I fl985). When the rotation of the white 
dwarf is close to break-up, numerical calculations showed 
that angular momentum will be transferred from the 
white dwarf into the accretion disk, and the rotation 
velocity will roughly keep constant (jPaczvnski 1 119911 : 
iPopham. fc Naravan I fl991) . Hence we set the specific 
angular momentum of the effective accreted matter to 
be r 

v < 0.9i>k, 



\J GM wd i?, 
, 



v > 0.9i>k 



(4) 



where v = u>R, and vk — \f GM wd /i? are the rotation 
velocity and the Keplerian velocity at the white dwarf's 
equator, respectively. In our calculations, the initial sur- 
face velocity at the white dwarf's equator is taken to 
be 10km s _1 , and the white dwarf radius changes with 

The losses of mass and orbital angular momentum play 
an important role in the mass transfer and orbital evo- 
lution of close binary systems. The mass loss rate of the 
binary system is M = (1 — X) a H,iC*He,i/i)M d , and we 
assume that this mass is ejected in the vicinity of the 
white dwarf in the form of isotropic winds or outflows, 
taking away the specific orbital angular momentum of 
the white dwarf. The rate of orbital angular momentum 
loss through mass loss is given by 



MM 



wcl 



(5) 



where J = (M d M wd /M)a 2 n, M = M d + M wd , and Q are 
the orbital angular momentum, total mass, and orbital 

angular velocity of the binary s ystem, respectivel y. 

With numerical simulations, lYoon fc Langer I (|2004f ) 
found that, accreting at a rate > 3 x 10 -7 M Q yr~ 1 , a 
white dwarf would rotate differentially, and there would 
not be the central carbon ignition even if its mass ex- 
ceed 1.4M . As a result of the differential rotation, 
a white dwarf with a super-Chandrasekhar mass may 
exist. Accordingly, we assume that a SN la occurs 
when M wd > 1.4M and M < 3 x lO _7 M yr" 1 , so 
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that there is no differential rotation to support the mas- 
sive white dwarf. In the case of M w d > 1.4M© and 
M > 3 X lO -7 M yr- 1 , we let the white dwarf increase 
mass up to M max - 

3. RESULTS 

We have calculated the evolution of white dwarf bi- 
naries adopting an upd ated version of the stel lar evolu- 
tion code developed by lEggleton I (11971L 11972ft (see als o 
lHan. Podsiadlowski fc Eggleton Ill994t iPols et atlll995ft. 
The s tellar OPAL opacities are from iRogers fc Iglesiasl 
(|1992ft . and lAlexander fc Ferguson I (|1994l) for alowtem- 
peratures. In the calculations we take the ratio of the 
mixing length to the pressure scale height to be 2.0. The 
evolutionary results are determined by three parameters 
of white dwarf binaries: the initial mass of the white 
dwarf M w d,i (taken to be 1.2M Q and 1.0M©), the initial 
mass of the donor star Md,i, and the initial orbital period 
of binary systems P rb,i- 

We first present an example of the evolutionary se- 
quences for a binary system with M w d,i = 1.2M©, M^i = 
2.5M©, and P rt>,i = 1-0 day, in which the donor star has 
a solar composition (Y = 0.28, Z = 0.02). In Figure 1 
we plot the evolution of the mass transfer rate and the 
donor star mass with time in the left panel, the evolu- 
tion of the orbital period and the white dwarf mass in 
the middle panel, and the evolution of rotation veloc- 
ity at the white dwarf's equator in the right panel. The 
donor star fills its Roche lobe when its age is ~ 3.84 x 10 s 
yr. Because of a relatively high initial mass ratio, the 
mass transfer occurs on a thermal timescale at a high 
rate of ~ 10~ 7 — 10~ 6 Af©yr _1 , but remains stable due to 
strong isotropic wind from the white dwarf. After ~ 1.5 
Myr mass transfer, the white dwarf grows to 1.75M©, 
and trigger a type la supernova when the mass trans- 
fer rate M declines to be 3 x 10 _7 M©yr _1 . The orbital 
period firstly decreases to 0.7 day, because material is 
transferred from the more massive donor star to the less 
massive white dwarf, and then increases when the white 
dwarf mass grows and exceeds the donor star mass. With 
gaining spin angular momentum from the accretion ma- 
terial, the rotation velocity of the white dwarf gradually 
increase, and near the break-up velocity. 

To investigate the distribution of the initial donor 
star mass and orbital period for the progenitor systems 
of super-Chandrasekhar mass SNe la, we have calcu- 
lated the evolutions of a large number of white dwarf 
binaries with different values of Md,i and P rb,i when 
Mvd,i = 1-2M©. Figure 2 summarizes our calculated 
results in the Mjj — P or b,i plane with Z = 0.02 (left 
panel) and 0.001 (right panel), respectively. In the case 
of Z = 0.02, the regions enclosed by the solid, dashed, 
and dotted curves denote the distribution areas of ini- 
tial white dwarf binaries with Msn > 1.7M©, > 1.6M©, 
and > 1.4M©, respectively. Beyond these areas, SNe 
la cannot occur due to either a low mass accumulation 
efficiency of the white dwarf or unstable mass transfer. 
In the case of Z = 0.001, the solid, dashed, and dot- 
ted curves correspond to the boundaries of distribution 
area for the progenitors with Msn > 1.6M©, > 1.5M©, 
and > 1.4M©, respectively, and we cannot find that a 
white dwarf mass exceed 1.7M©. When the metallicity 
Z changes from 0.02 to 0.001, the occupied region by 
the progenitors of SNe la has a tendency to move down- 



ward in the Md,; — P rb,i diagram, i.e., the progenitors 
with higher Z t end to have a more massive donor star, in 
agreement wi th lMeng. Chen fc Han 1 (2009). However, it 
is difficult to produce a super-Chandrasekhar mass SN la 
when the donor mass is low as 1.0M©. Figure 3 presents 
the progenitor distribution of SNe la in the Md,i — P rb,i 
plane when M w d,i = 1.0M©. For Z = 0.02, the maxi- 
mum explosion mass of white dwarfs is always less than 
1.5M©. 

Figure 4 shows the distribution of A/d.f and P rb,f at 
the moment of SN explosion. The open triangles, open 
circles, and solid stars denote systems with 1.6M© > 
Msn > 1-4M©, 1.7M© > M SN > 1.6M©, and M SN > 
1.7M© in the left panel (Z = 0.02), and 1.5M© > M SN > 
1.4M©, 1.6M© > M SN > 1.5M©, and M SN > 1.6M© in 
the right panel (Z = 0.001), respectively. For Z = 0.02 
and Msn > 1.7M©, the companion masses ~ 1.3— 1.8M© 
and orbital periods ~ 0.4 — 2.0 d, while in the case of 
Z = 0.001 and Msn > 1.6M©, the companion masses 
~ 0.5 — 1.5M© and orbital periods ~ 0.4 — 6.0 d. In 
Figure 5, we plot the distribution of the companion stars 
at the moment of SNe la in the H-R diagram. They have 
luminosities in the range of ~ 1.0 — 10L© and effective 
temperatures in the range of ~ 5000 — 6000 K. These 
properties may be compared with and testified by future 
optical observations of SN la remnants. 

Certainly our approach in estimating the mass accumu- 
lation efficiency on rotating white dwarfs is very rough 
and simplified. We only consider the effect of centrifu- 
gal force and neglect the thermal and chemical evolution 
on the surface of rapid rotating white dwarfs, which is 
very complicated. In our calculations we assumed a rigid 
rotation instead of a differential rotation, which only in- 
crease slightly the Ch andrasekhar limit to be 1.48M© 
(|Yoon fc Langer 1 12004| ). Additionally, the white dwarf 
will be deformed during spin-up, not obeying spherical 
symmetry as we assumed. Nevertheless, Our results show 
that the mass accumulation efficiency of white dwarfs de- 
crease wit h rotation. This is at le ast qualitatively consis- 
tent with iPiersanti et al. I (|2003T ) , in which the authors 
compared rigidly rotating white dwarfs with non-rotating 
objects for CO-accretion at various accretion rate, and 
found that an increase in the angular velocity of the 
white dwarf always leads to a decrease in the value of 
the accretion rate below which central carbon ignition 
will occur. The reason is that including rotation cause 
the accreting star to be both less dense and cooler, re- 
sulting in a larger thermal diffusion timescale, while the 
effect of compression induced by the accretion process is 
only slightly modified. This implies that our treatment 
might not be far from real situations. 

The influence of metallicities on the mass accumula- 
tion efficiency can be found in Figure 6, which shows the 
fraction of mass growth in the transferred mass 

M w d,f — M w d,i , . 



against the initial orbital period P rb,s when the initial 
mass of the donor star Md.i = 2.5M©. It is clearly seen 
that, for the same donor star, white dwarfs accompanied 
by a Population I star have a higher mass growth rate. 

4. DISCUSSION AND SUMMARY 
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Recently discovered overluminous SNe la sug- 
gested that they might have ori ginated from super- 
Chandrasekhar mass white dwarfs (IHowell et al. I l2006t 
iHicken et all l2007t lYuan et al. 1 |2007|). Based on 
the single degenerate model and the suggestion 
that massive white dwarfs might be supported by 
rapid differential rotation when accretion rate M > 
3 x l Q-^^vr- 1 (lDurisen1 ll97a IS napiro fe Teukolskvl 
l!983HYoon fe Langer II2005I ). we have performed numer- 
ical calculations of the evolution of white dwarf binaries 
to investigate the properties of the progenitor of super- 
Chandrasekhar mass SNe la. The results can be summa- 
rized as follows. 

(1) For white dwarfs with an initial mass of < I.OMq, 
the explosion masses of SNe la are nearly uniform, i.e. 
super-Chandrasekhar mass SNe la are difficult to pro- 
duce. This is consistent with the rareness of super- 
Chandrasekhar mass SNe la in observations. 

(2) When M w a,i — 1-2M©, depending on the evolution- 
ary paths the masses (-Msn) of exploding white dwarfs 
range from 1.4Af to 1.76A/q. A considerable fraction 
of SNe la are of super-Chandrasekhar mass, suggest- 
ing a diversity in the brightness of SNe la. However, 
in most cases the final masses of the white dwarfs are 
not significantly exceeding 1.4 Mq, and it is very diffi- 
cult to produce a SN la with Msn > 1-8 Mq. Thus 
our model cannot reproduce an ov erluminous SN la wi th 
M SN > 2.0Mq like SNLS-03D3bb (IHowell et al. Il2006h 2 . 

(3) Progenitors of super-Chandrasekhar mass (M SN > 
1.6M Q ) SNe la can be constrained to be white dwarf 
(with initial mass of l.2M@) binaries with an initial 
donor star mass of Md,i ~ 2.2 — 3.3M© and an initial 
orbital period of P orb ,i ~ 0.5 - 4.0 d when Z = 0.02, 
and M d>i ~ 1.7 - 2.7M Q and P orb ,i ~ 0.5 - 3.5 d 
when Z = 0.001. Interestingly, it is not those with 
most massive donors that are more likely to evolve to 
super-Chandrasekhar mass SNe la, since the mass trans- 
fer rates in these systems usually decline rapidly below 
3 x 10 -7 Af©yr -1 , so that the white dwarf does not have 
sufficient time to accrete. The time delay between the 
formation of the progenitor systems and the explosions of 
SNe la is < 1 Gyr, in agreement with the suggestion that 
super-Chandrasekhar mass SNe la should be more likely 
to ex ist in a young stellar population ([Howell et al. I 
2006), but they seem not to belong the youngest pop- 

2 It is controversial if the progenitor of SN 2003fg is a super- 
Chandrasekhar mass white dwarf. Based on the low velocity and 
short time-scale seen in SN 2003fg, Macda & Iwamoto (2009) sug- 



ulation. Especially one would not expect overluminous 
SNe la in early type galaxies. 

(4) Metallicities have important influence on the pro- 
duction of super-Chandrasekhar mass SNe la. Systems 
with Population II donor stars are less likely to be the 
progenitor of SN la with a super-Chandrasekhar mass 
of > 1.7M©, suggesting that SNe la at relatively high 
rcdshifts might be more homogeneous than nearby ones. 

Perhaps the biggest issue in this work is how super- 
Chandrasekhar mass white dwarfs can be formed and 
related to overluminous SNe la. The nature of the over- 
luminous SNe la is not yet well understood. For example, 
Subaru and Keck optical spectroscopic and photometric 
observations of SN la 2006gz show that, the late-time 
behavior of this SN is distinctly different from that of 
normal SNe la, but the peculiar features found at late 
times (the SN is faint and it lacks [Fe II] and [Fe III] 
emiss ion) are not readily connected to a large amount of 
56 Ni (jMaeda et al. Il2009h . Even if super-Chandrasekhar 
mass white dwarfs do exist in the universe, their for- 
mation path is poorly known. In principle, both merge 
of a close binary system of two massive white dwarfs 
and accretion onto a normal white dwarf may form a 
rapidly spinning white dwarf that considerably exceeds 
the Chandrasekhar limit, provided that they rotate dif- 
ferentially. In the latter case, under what circumstances 
can a differentially rotating white dwarf be sustained has 
not been well addressed. We have adopted the simple 
criterion of the mass accretion rate, but the real situa- 
tion should be sensitively dependent on the efficiency of 
mass and angular momentum gain and loss during mass 
transfer, which may be also related to the magnetic field 
strength of the white dwarfs. 
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Fig. 1. — Evolution of a white dwarf binary with Af wt j ; = 1.2Mq, i = 2.5Mq and P or b i = l.Oday. The solid and dashed curves 
denote the evolutionary tracks of the mass transfer rate and the donor star mass in the left panel, the orbital period and the white dwarf 
mass in the middle panel, and the rotation velocity and break-up velocity at the white dwarf's equator, respectively. 
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Fig. 2. — Distribution of the initial orbital periods P or b.i and the initial donor star masses Md,! of the progenitor systems of supcr- 
Chandrasekhar mass SNe la when M wt j_i = 1.2Mq. The left and right panels are for Z = 0.02 and 0.001, respectively. 




Fig. 3. — Distribution of the initial orbital periods P or b.i an d the initial donor star masses M^^ of the progenitor systems of SNe la when 
A^wd.i = 1-OMq. The left and right panels arc for Z = 0.02 and 0.001, respectively. 
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Fig. 4. — The donor star masses Md,f vs. the orbital periods P rb.f when SN la explosions occur. In the left panel (Z = 0.02 ), the 
open triangles, open circles, and solid stars denote systems with 1.6Mq > Msn > IAMq, 1.7Mq > MgN > 1.6Mq, and Msn > 1.7 Mq, 
respectively. In the right panel (Z = 0.001), the open triangles, open circles, and solid stars denote systems with 1.5Mq > Mg^ > 1.4MQ, 
1.6M > M SN > 1.5M , and M SN > 1.6M , respectively. 
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Fig. 5. — Distribution of the donor stars in the H-R diagram when SN la explosions occurs. Symbols are the same as in Figure 4. 
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Fig. 6. — The parameters rj vs. the initial orbital periods P rb,i when the initial masses of donor stars Mj.; = 2.5Mq. The filled circles 
and open circles correspond to the calculated results for Z = 0.02 and 0.001, respectively. 



